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Abstract 

Background: Air travel has expanded at an unprecedented rate and continues to do so. Its effects have been seen 
on malaria in rates of imported cases, local outbreaks in non-endemic areas and the global spread of drug 
resistance. With elimination and global eradication back on the agenda, changing levels and compositions of 
imported malaria in malaria-free countries, and the threat of artemisinin resistance spreading from Southeast Asia, 
there is a need to better understand how the modern flow of air passengers connects each Plasmodium 
falciparum- and Plasmodium wvw-endemic region to the rest of the world. 

Methods: Recently constructed global P. falciparum and P.vivax malaria risk maps, along with data on flight 
schedules and modelled passenger flows across the air network, were combined to describe and quantify global 
malaria connectivity through air travel. Network analysis approaches were then utilized to describe and quantify the 
patterns that exist in passenger flows weighted by malaria prevalence. Finally, the connectivity within and to the 
Southeast Asia region where the threat of imported artemisinin resistance arising is highest, was examined to 
highlight risk routes for its spread. 

Results: The analyses demonstrate the substantial connectivity that now exists between and from malaria-endemic 
regions through air travel. While the air network provides connections to previously isolated malarious regions, it is 
clear that great variations exist, with significant regional communities of airports connected by higher rates of flow 
standing out. The structures of these communities are often not geographically coherent, with historical, economic 
and cultural ties evident, and variations between P. falciparum and P. vivax clear. Moreover, results highlight how 
well connected the malaria-endemic areas of Africa are now to Southeast Asia, illustrating the many possible routes 
that artemisinin-resistant strains could take. 

Discussion: The continuing growth in air travel is playing an important role in the global epidemiology of malaria, 
with the endemic world becoming increasingly connected to both malaria-free areas and other endemic regions. 
The research presented here provides an initial effort to quantify and analyse the connectivity that exists across the 
malaria-endemic world through air travel, and provide a basic assessment of the risks it results in for movement of 
infections. 



Background 

The worldwide air travel network has expanded at an ex- 
ceptional rate over the past century. International pas- 
senger numbers are projected to rise from 1.11 billion in 
2011 to 1.45 billion by 2016, with an annual growth rate 
of 5.3% [1]. Today, there are 35,000 direct scheduled 
routes on the air travel network, with 865 new routes 
established in 2011 [2]. Malaria-endemic areas are more 
connected to the rest of the world than at any time in 
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history, with the disease able to travel at speeds of 600 
miles per hour within infected passengers. The growth 
of the air travel network results in substantial concerns 
and challenges to the global health system, with a need 
to place more emphasis on evidence-driven surveillance 
and reporting that incorporates spatial and network in- 
formation [3-6]. 

Rising rates of travel between malaria-free and -endemic 
countries have led to general patterns of increased rates of 
imported malaria over recent decades [7-10]. Due to infre- 
quent encounters [3,9], imported cases can challenge 
health systems in non-endemic countries, with difficulties 
in diagnosis [10], misdiagnosis and delays in treatment 
[11,12], as well as significant treatment expenses [13]. 
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Further, flights may bring infected vectors, resulting in "air- 
port malaria", where patients who do not have a foreign 
travel history become infected through being bitten in the 
vicinity of international airports [14-17]. Patterns in 
imported cases and airport malaria have been shown to be 
related to a combination of the numbers of travellers and 
the malaria risk at the destination [3,16], and these relation- 
ships will continue to evolve as new routes become 
established. 

The flow of people via air travel between endemic areas 
may increase the risks of re-emergence or resurgence [18] 
in previously malaria free or low transmission areas [19]. 
The autochthonous malaria outbreaks in Virginia in 2002 
[20], Florida in 2003 [21] and Greece in 2011 [22], for ex- 
ample, demonstrate the continued risks of local outbreaks 
following reintroduction through air travel, though such 
occurrences are rare [23]. Further, the examples of malaria 
resurgence in island nations, such as Sri Lanka [24], 
Mauritius [25] and Madagascar [26], after control mea- 
sures were relaxed reinforce the importance of vigilance 
and robust surveillance in terms of human movement in 
pre and post-elimination periods [18]. Identifying the risks 
of malaria movement through the air travel network can 
provide an evidence base through which public health 
practitioners and strategic planners can be informed about 
potential malaria influxes and their origins [3,27]. 

Meanwhile, growing concerns have been raised about 
the possible spread of artemisinin resistance from the 
Greater Mekong subregion in Southeast Asia to other 
endemic areas. Recent research has highlighted increasing 
numbers of patients showing slow parasite clearance rates 
following treatment with artemisinin-based drugs in 
the Cambodia-Thailand border and Thailand-Myanmar 
border regions [28-30]. Tremendous health and socio- 
economic costs occurred when chloroquine-resistant par- 
asites arrived in sub-Saharan Africa from Southeast Asia 
and spread across the continent [31,32]. Similarly, 
sulphadoxine and pyrimethamine resistance emerged in 
Asia and spread to Africa [33,34]. The WHO reports that 
there is already "at least one study with a high treatment 
failure rate (>10%) reported from six of the 23 African 
countries that have adopted artesunate-amodiaquine com- 
pound" [35], and fear remains over the spread of artemisi- 
nin resistance from Southeast Asia to Africa, that could 
undermine current control and elimination efforts, with 
no alternative drugs coming in the foreseeable future. 

Rates of imported malaria, risks of resurgence and the 
spread of drug resistance are all today influenced by how 
the global air travel network connects up the malaria- 
endemic regions of the world, and the numbers of passen- 
gers moving along it. Here, recently constructed global 
Plasmodium falciparum and Plasmodium vivax malaria 
prevalence maps are combined with data on modelled pas- 
senger flows across the air network, to describe and 



quantify global malaria connectivity through air travel in 
2010. Weighted network analysis statistics are derived to 
examine: (i) which regions show greatest connectivity to 
P. falciparum and P. vivax malaria-endemic zones; (ii) where 
the largest estimated passenger flows from endemic areas 
occur; (iii) which regions form communities; whereby mal- 
aria infection flows within them are likely to be larger than 
between communities, and finally, (iv) where the threat of 
imported artemisinin resistance is highest via air traffic, 
and the possible risk routes for the spread of resistance 
within and from Southeast Asia. 

Methods 

Airport locations, flight routes and passenger flow matrix 

Information on the longitude, latitude, city name and air- 
port code for a total of 1,449 airports which serve cities 
with more than 100,000 people, and a modelled actual' 
traffic flow (i.e. number of passengers travelling between 
each location and every other, irrespective of stopovers) 
connectivity list with 644,406 routes amongst these air- 
ports were obtained [6,36,37]. This list documented, for 
each origin and final travel destination, the estimated 
number of passengers taking this route [36] regarding the 
hub-and-spoke structure of the air travel network [38]. A 
connectivity matrix was then created from the connectiv- 
ity list, quantifying the volumes and the directionalities of 
the passenger flows between two airports. Within this pas- 
senger flow matrix, 23,785 origin-destination pairs were 
connected by direct flights between two airports, 291,745 
pairs were connected by routes involving one-stop and 
328,876 pairs required two-stops to connect. The travel 
volumes on the routes were modelled based primarily on 
publicly available datasets under a generalized linear 
model framework. Full model details are provided in 
Huang et al [36], but in brief, to construct the matrix, 
topological characteristics of the air travel network, city 
population, and local area GDP, amongst others, were uti- 
lized as covariates. Actual travel volumes for training and 
validation were extracted and assembled from various 
transportation organizations in the USA, Canada and the 
European Union. A log linear model controlling for ran- 
dom effects on origin, destination and the airport hier- 
archy was then built to predict passenger flows on the 
network. The model outperformed existing air travel pas- 
senger flow models in terms of prediction accuracy [36]. 

Malaria distribution 

Global P. falciparum and P. vivax prevalence maps were 
obtained from the Malaria Atlas Project [39] and the 
methods behind their construction are presented in 
Gething et al. [40,41]. In brief, 22,212 community preva- 
lence surveys were used in combination with model- 
based geostatistical methods to map the prevalence of 
P. falciparum globally in 2010 within limits of transmission 
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defined by annual parasite incidence and satellite covari- 
ate data. Similarly, 9,970 geocoded P. vivax parasite rate 
(P. vivax PR) surveys collected between 1985 and 2010 
were utilized in a spatiotemporal Bayesian model-based 
geostatistical approach to map endemicity, under the re- 
strictions of a mask of the stable/unstable endemicity 
[41] and information on the prevalence of the Duffy 
blood group [42]. Distributions of Plasmodium ovale, 
Plasmodium malariae or Plasmodium knowlesi are not 
included here, since similar datasets on their distribu- 
tions do not yet exist. Also, the seasonal climatic con- 
straints that affect the transmission of P.falciparum and 
P. vivax are not included here, but models of each [43] 
will be included in future work. 

Weighted network analysis and community detection 

Malaria prevalence can vary greatly in the region around 
airports and the cities they serve, and travellers taking 
flights from a specific airport may reside many 
kilometres from the airports in higher transmission areas 
than found in the vicinity of the airport. Thus, simply 
assigning the predicted prevalence from the malaria 
maps at the location of each airport could underestimate 
the risk and rate of infection exportation at the airport 
in question and underrepresent its contribution to global 
malaria connectivity. Therefore, following Huang et al 
[6], local accessibility to each airport was considered by 
assuming that passengers would travel less than 50 km 
with a travel time less than two hours to access an 
airport to take a flight. Under this assumption, the 
P. falciparum and the P. vivax malaria prevalence assigned 
to an airport were obtained as the maximum prevalence 
from the malaria maps within a mask of 50 km and two- 
hour travel time (Additional file 1), in which the mask 
was generated using a global travel time map [44]. This 
choice of the maximum value was created for a worst 
case scenario' assessment, but for most of the analyses 
conducted here, relative differences between airports 
were assessed or presented, thus making the choice of, 
for example, maximum, minimum or mean, irrelevant in 
most cases. Likewise, an indicator that defines whether 
an airport is located in the stable/unstable endemic zone 
was created according to the same mask, in which the 
indicator defines whether the majority area of the mask 
is located in the stable/unstable zone. Additional file 1 
shows this travel time/distance mask with the global 
travel time map. 

The above approach ensured that each airport had an 
assignment of a P. falciparum and P. vivax prevalence 
rate (or unstable/malaria free), which could then be used 
as a weighting applied to the passenger flow estimates to 
derive relative malaria flow' indices for each route, 
that could be compared to other routes across the 
global network to analyse malaria connectivity. Thus, 



P. falciparum and P. vivax flows were calculated on each 
route (either direct, one-stop or two-stops) as origin 
prevalence * estimated passenger volume, to produce 
P. falciparum and P. vivax malaria networks. 

A group of weighted centrality analyses and network 
community partition analyses were performed on the 
malaria networks to quantify features of global malaria 
connectivity. First, the in-strength and out-strength of 
each connection was calculated as the sum of incoming 
and outgoing malaria flow via all possible connections 
(direct flight, one-stop or two-stops) as follows: 

N 

s. = > a.,w„ (1) 

i Z ✓ i) ij v 7 

7=1 

in which a t j is the airport adjacency matrix (in a binary 
form) and is the weighted malaria flow. This metric 
estimates the total weight of malaria flows that airports 
send and receive. 

Following this, weighted "betweenness" analyses were 
performed on the malaria flow matrices. Betweenness cen- 
trality measures the number of shortest paths going 
through a specific vertex [45]. In a weighted network, be- 
tweenness centrality is a useful local measure of the load 
placed on the given node in the network as well as the 
node s importance to the network other than just connect- 
ivity [46]. It is often used in transport network analysis to 
provide an approximation of the traffic handled by the 
vertices [47]. Thus, here it provides an indication of the 
status of each airport as a 'malaria hub' through its im- 
portance in the global flow of malaria infections via air 
travel - i.e., a measure of how many infections likely pass 
through each airport each year, relative to other airports, 
and how likely an airport would route and spread malaria 
infections onward. The betweenness centrality is calcu- 
lated as: 

Cb(v) = E^ W 

in which S st is the total number of shortest paths from 
node s to node t and S st (v) is the number of those paths 
that pass through v. Note that on the weighted P. falcip- 
arum/P. vivax networks, the distance between the two 
nodes s and t is defined by the sum of P. falciparum flows 
or P. vivax flows as the edge weight on this path under a 
classical Dijkstra shortest-path calculation framework [48]. 
A normalized betweenness was used as 

CBnormiy) = v^vT^ \ ~ ( 3 ) 

^C B (v)/n 

where n was the number of nodes (airports) in the air 
travel network [47] . 
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Communities in a network reflect the partition of 
nodes that are densely connected and separated from 
the other nodes in the network, thus these nodes "prob- 
ably share common properties and/or play similar roles 
within the graph" [49]. By mapping communities on the 
malaria networks defined here, groups of airports that 
show strong links in terms of likely movements of infec- 
tions were identified. This potentially has utility in terms 
of providing evidence upon which regional surveillance 
strategies can be designed [50,51]. Newman and Girvan 
[52] define a modularity score which measures the qual- 
ity of network partitions as: 



2m 4 



2m 



S{c h cj) 



(4) 



in which, Wy represents the weight of the edge between i 
and j (here these are the P. falciparum and P. vivax flow 
matrices), k t = X jW^ is the sum of the weights of the 
connections attached to airport i, c t is the community to 
which airport i is assigned; £fec ; ) is 1 if c { = Cp otherwise 

A multilevel algorithm for community detection [53] 
was implemented. This method utilizes an iterative ap- 
proach that merges communities to maximize the 
modularity score: Firstly, modularity is optimized by 
allowing only local changes of communities; secondly, 
the established communities are combined together to 
construct a new network. These two passes are repeated 
iteratively until no increase of modularity is possible. 
The number of communities returned by this algorithm 
yields the maximum modularity score. 

Afterward, a simple Wilcoxon rank-sum test [54] was 
performed on the differences between "internal" and 
"external" degrees of a community in order to test 
whether the establishment of communities was signifi- 
cant. Air connections were defined within a community 
as "internal" and the connections connecting the air- 
ports of a community with the rest of the network as 
"external". The null hypothesis of this test was that there 
was no difference between the number of internal and 
external routes incident to an airport of the community. 

Results 

The results of the global malaria connectivity analyses 
are presented in two sections: (i) analyses focussed on 
the connection of endemic malaria regions to each other 
and to malaria-free areas, that has particular relevance 
to imported malaria and malaria resurgence and re- 
emergence; and, (ii) analyses examining the connections 
between Southeast Asia and the rest of the malaria- 
endemic world, which are relevant to the spread of arte- 
misinin resistance. 



Connectivity within endemic areas and to non-endemic 
areas 

Figure 1 shows the results of regional community struc- 
ture analyses based on traffic flow data overlaid on the 
P. falciparum/P. vivax endemicity and stable/unstable 
transmission limits maps. The Wilcoxon test results 
show that the internal degrees for the airports within all 
communities are significantly different from the external 
degrees, with p values of <0.01, thus the community par- 
titions shown are significant. The maps highlight those 
countries that form communities linked by high levels of 
traffic scaled by P. falciparum/P. vivax prevalence at 
their origin endemic area. Additional file 2 describes 
similar analyses based solely on the travel network data 
from Huang et al [36]. The communities detected re- 
flect the architecture of the air network, and how this 
relates to malaria endemicity around the world. Geo- 
graphical contiguity is clearly evident, as traffic levels on 
shorter distance routes are generally higher than on lon- 
ger distance routes, but interesting patterns relating to 
historical ties emerge. For instance, for P. falciparum, 
London forms part of the Nigeria community, but Paris 
shows stronger ties to the remainder of sub-Saharan 
Africa. These connections are often reflected in im- 
ported malaria statistics, with Nigeria being the main 
source of P. falciparum cases seen in the UK, but for 
France, the French-speaking African countries are the 
main origin. Similarly, UK airports also form part of the 
India/Bangladesh community, where historical ties exist, 
resulting in significant travel between the two regions, 
and consequent P. falciparum and P. vivax malaria im- 
portation to the UK. Ties also exist between the western 
USA and East Asia, which form a single P. falciparum 
community (Figure 1A). Additional file 3 shows a com- 
munity detection analysis for airports with direct con- 
nections, one-transfer connections and two-transfer 
connections from malaria-endemic areas. 

To examine directional and net potential movements 
of people and parasites between airports in different 
countries, the international route weightings were 
summed up to identify possible "source" and "sink" air- 
ports of malaria infections (Additional file 4). Here, the 
weights of all possible incoming flows for airports in the 
non-endemic areas, and the weights of all possible out- 
going flows from airports in endemic areas were 
summed up to define "vertex strengths" of importation 
and exportation (note that only the routes connecting 
two different countries regardless of the domestic routes 
were considered). In this table, airports in the Far East 
and Middle Asia such as Singapore, Hong Kong, Dubai, 
and Sharjah display the highest importation values (note 
that Singapore ranked the first in both categories). Un- 
surprisingly, major air hubs in Europe (such as airports 
in London, Paris and Frankfurt) also showed high 
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(See figure on previous page.) 

Figure 1 Spatial distribution of Plasmodium falciparum/Plasmodium vivax network communities overlaid on Plasmodium falciparum/ 
Plasmodium vivax prevalence maps. A) P. folciporum multilevel membership. B) P. vivox multilevel membership. These two maps show only 
airports that have direct connections from endemic to non-endemic areas, though full origin-destination flow estimates were used in 
calculations. The inset maps present close-up views of the USA and western Europe. Airports with the same community membership (indicated 
by the same colour) display stronger links in terms of likely movements of infections between them than to airports in other communities. Note 
that in the P. vivax map, communities with less than ten airports are not shown. 



potential incoming P. falciparum flows. Miami is the 
only airport in the USA on the importation flow top ten 
lists, with its strong connections to Central and South 
America. In terms of exportation, the largest airports by 
traffic capacity and connections to the rest of the 
malaria-endemic world were highlighted (Additional file 4). 
Mumbai was ranked first as the largest exporter of 
P. falciparum and P. vivax flows, suggesting that it likely 
acts as an important portal for spreading malaria to the rest 
of the world. 

The betweenness centrality metric was utilized to in- 
spect the connectivity from endemic areas. As the be- 
tweenness metric C B (v) is defined as the number of 
shortest paths connecting any two airports that involve a 
transfer at airport v, high centrality airports in endemic 
areas provide hubs for people originating at less- 
accessible airports in remote places to reach the rest of 
the world. Thus the betweenness centrality presents an 
initial quantification for the likelihood that an airport 
routes or spreads malaria flow as a 'malaria hub'. Table 1 
shows the top ten highest betweenness centrality air- 
ports for transferring P. falciparum flow and P. vivax 
flow elsewhere. For the P. falciparum flow, international 
airports in Africa play important roles as hubs for rout- 
ing infections. Some airports are observed to have small 
degrees (low numbers of connecting routes) and large 
centrality (importance as a hub), which can be consid- 
ered as an abnormality [47]. These airports connect less 
accessible and connected airports in endemic areas to 
other airports in the world. For P. vivax flow, Asian 
international hubs play more important roles. Of interest 
is Phoenix airport, which ranked the sixth in terms of 
P. vivax centrality, suggesting that it plays an important 
role as a gateway in linking P. Wv<2#-endemic areas to 
the USA. Additional file 5 presents the spatial distribu- 
tion of betweenness centrality scores for airports, 
weighted by P. falciparum or P. vivax flows. 

To further investigate the effects of flows from endemic 
zones, Additional file 6A and B shows the sums of inter- 
national incoming risk flows for all the airports in those 36 
countries that have national policies for malaria elimin- 
ation, and are closest to eliminating the disease [55]. Im- 
portation of infections threatens the success of elimination 
programmes [19] and while air travel may not be the 
highest risk source for these introductions for most of 
these countries, it remains a potentially important source 



of incoming infections. From these two maps, it can be 
seen that China and countries in Middle Asia are subjected 
to the greatest pressure of incoming flows, relative to other 
elimination countries, due to their larger incoming traffic 
volumes from endemic regions elsewhere around the 
world. To sum up, detailed analyses on airport connectivity 
are provided in Additional files 2, 3, 4, 5 and 6. 

Connectivity to Southeast Asia 

Figure 2 maps out the passenger flows scaled by origin 
prevalence for P. falciparum and P. vivax from the 
Greater Mekong subregion. Significant amounts of flow 
exchange within Southeast Asia can be seen in the 
close-up subsets. For both P. falciparum and P. vivax it 
can be seen that the connectivity, through numbers of 
travellers, to Latin American endemic regions is weak, 
but that much stronger connections to sub-Saharan 
Africa and the Indian subcontinent exist. Increasing con- 
nections through trade and labour markets between Asia 
and Africa over the past decade is exemplified here in 
the strong connections between the Southeast Asian re- 
gion and all of sub-Saharan Africa's major airport hubs. 
Additional file 7 presents the top ten risk routes spread- 
ing drug resistance of P. falciparum and P. vivax from 
the Greater Mekong subregion to non-Asian destina- 
tions, with estimated P. falciparum /P. vivax flow and 
the number of stops needed to travel from the origin city 
to the destination city shown. 

Discussion 

The continuing growth in air travel is playing an import- 
ant role in the global epidemiology of malaria. Flight 
routes now connect previously isolated malaria-endemic 
regions to the rest of the world, and travellers on these 
routes can carry infections to the opposite side of the 
world in less than 24 hours. While many endemic areas 
still remain relatively isolated, the malaria-endemic 
world is becoming increasingly connected to both 
malaria-free areas and other endemic regions. The im- 
pacts of this can be seen in imported cases, vector inva- 
sions and the spread of drug-resistant parasite strains. 
Here a spatial network analysis approach were presented 
to demonstrate the connectivity that exists across the 
malaria endemic world through air travel, and provide 
quantitative indicators of the risks it results in for mal- 
aria movement. 
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Table 1 Top ten Plasmodium falciparum/Plasmodium vivax betweenness centrality airports with their degrees in a 
sub-network that only contains direct links from airports in Plasmodium falciparum- or Plasmodium wVax-endemic areas 

Top ten P. falciparum centrality airports 



Airport 


City 


Country 


Normalized betweenness centrality 


Degree 


NBO 


Nairobi 


Kenya 


47.35 


80 


MBA 


Mombasa 


Kenya 


32.44 


27 


JRO 


Kilimanjaro 


Tanzania 


32.39 


14 


BOM 


Mumbai 


India 


30.41 


104 


ADD 


Addis Ababa 


Ethiopia 


28.21 


64 


DEL 


Delhi 


India 


23.16 


111 


JIB 


Djibouti 


Djibouti 


19.77 


15 


ADE 


Aden 


Yemen 


18.63 


15 


MGQ 


Mogadishu 


Somalia 


14.45 


8 


HRE 


Harare 


Zimbabwe 


14.35 


20 


Top ten P. vivax centrality airports 


Airport 


City 


Country 


Normalized betweenness centrality 


Degree 


BKK 


Bangkok 


Thailand 


96.43 


146 


ICN 


Seoul 


South Korea 


78.12 


150 


DEL 


Delhi 


India 


59.55 


133 


BOM 


Mumbai 


India 


34.17 


116 


KMG 


Kunming 


China 


30.79 


90 


PHX 


Phoenix 


USA 


28.63 


91 


DPS 


Denpasar Bali 


Indonesia 


27.94 


34 


SJO 


San Jose 


Costa Rica 


27.72 


37 


DOH 


Doha 


Qatar 


25.91 


100 


TAS 


Tashkent 


Uzbekistan 


25.85 


69 



The betweenness centrality scores show how many shortest paths go through an airport, thus they highlight the potential that airports might route infection 
flows, acting as 'malaria hubs'. The degree measures how many routes are linked to an airport. 



Results highlight the substantial connectivity that now 
exists between and from malaria-endemic regions 
through air travel While the air network provides con- 
nections to previously isolated malarious regions, it is 
clear that great variations exist, with significant regional 
communities of airports connected by high rates of 
prevalence-scaled flow standing out (Figure 1). The 
structures of these networks are often not geographically 
coherent, with historical, economic and cultural ties evi- 
dent. As new routes continue to be established, these 
communities will likely change, with new popular travel 
routes, such as those between China and Africa [3] likely 
altering global malaria flow routes, and new destinations 
that might encounter increased risks of imported mal- 
aria will emerge (Additional file 4 and Table 1). These 
community maps (Figure 1) and lists of cities by likely 
import/export of infections (Additional file 4) and hubs 
for infection flow (Table 1) provide a quantitative picture 
of how malaria infections are likely moving globally 
through air travel, and information from which global 
surveillance strategy design can draw upon. Additional 
file 4 and Table 1 highlight that certain airports provide 



significant hubs and gateways for the movement of in- 
fections and their entry into countries, and that these 
are widely distributed across the world. Their role in 
providing important nodes as both significant through- 
flow of infections in the network, and entry and exit 
gateways for cases to/from regions means that they po- 
tentially represent valuable sentinel sites for focussed 
surveillance. Finally, Figure 2 provides a stark reminder 
of how well connected the malaria-endemic areas of 
Africa are now to Southeast Asia, illustrating the many 
possible routes that artemisinin-resistant strains could 
take. These routes can provide a first-step quantification 
to support the global plan against artemisinin resistance 
containment [35] and design of surveillance systems 
[56], and should be refined with information on the lo- 
cations of resistance found. Such data could also inform 
decisions on where and how to limit the risk of spread, 
for example by pre-travel or arrival screening and 
treatment. 

A range of limitations and uncertainties exist in the 
analyses presented here. In terms of the quantification of 
malaria transmission, the use of static maps of annual 
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B 

Figure 2 Estimated relative Plasmodium falciparum/Plasmodium vivax flows originating from the Great Mekong subregion overlaid on 

Plasmodium falciparum/Plasmodium vivax prevalence maps. A) P. falciparum flows originating from the Great Mekong subregion. B) P. vivax flows 

originating from the Great Mekong subregion. The flows include estimated passenger numbers, including direct, one-transfer and two-transfer flight 

routes. The inset maps show close-up views for airports in the Greater Mekong subregion. 
) 



average prevalence [40,41] neglects the seasonality in 
transmission that is common to many areas, and also the 
substantial changes in transmission intensity seen in a var- 
iety of locations in recent years [57]. In most parts of the 
world, the densities of Anopheles mosquitoes change sea- 
sonally, thus impacting the receptivity of these areas to 
malaria flows incoming through air travel. While data on 
changes in Anopheles densities globally are not available, 
temperature-driven models of malaria transmission sui- 
tability [43] could be integrated in future work to better 



account for this and the seasonally varying topologies of 
the global malaria connectivity network studied. The 
demographic and behavioural differences between passen- 
gers are not accounted for here. Those taking regular air 
travel are often richer [58], and less likely to be infected, 
while those that are actually infected and showing symp- 
toms may be less likely to travel. Hence, the air travel 
passenger dataset used here clearly contains some biases 
when addressing malaria risks. Further, only parasite 
prevalence was used as a malaria metric, and while this 
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may be an adequate measure of population prevalence at ori- 
gin locations, it is not so appropriate for assessing the risk of 
infection acquisition for naive travellers, and entomological- 
based indices are likely more appropriate here, as used in 
more local studies [27,51,59]. Finally, the examination of rela- 
tive artemisinin resistance spread risk focuses simply on all 
travel from four countries, and thus does not account for 
any heterogeneity in resistance in the region. 

Uncertainties and limitations relating to the travel data 
used also exist. The modelled passenger flows represent 
just a 2010 snapshot, and thus routes and changes since 
then are not captured, while inherent uncertainty due to 
the modelling process also exists [36]. Moreover, the 
types of traveller and their activities during travel and 
their residential location are unknown, each of which 
contributes to differing malaria infection risks. Finally, 
overland and shipping travel flows are not considered 
here, which also contribute to local, regional and global 
malaria connectivity and flows. 

This work forms the basis for future analyses on 
imported malaria, elimination feasibility and the risks 
and potential routes of artemisinin resistance spread. 
Rates and routes of imported malaria have been shown 
to be significantly related to a combination of numbers 
of travellers to/from endemic destinations and the 
prevalence of malaria there [3]. The potential thus exists 
to construct a model based on global malaria prevalence 
[40,41], the local spatial interaction and accessibility to 
an airport within a region [60], transmission models for 
attack rate estimation [27], and traveller flow data [36], 
that can be used to forecast imported malaria rates, vali- 
dated with imported malaria data reported by health fa- 
cilities/organizations. 

As nations make progress towards elimination [55], 
the importance of human movement and imported cases 
increases. This work contributes to an on-going initia- 
tive, the human mobility mapping project [61], aimed at 
better modelling human and disease mobility, and will 
form one aspect of continued multimodal assessments 
of malaria movements [19,27,50,51] and assessment of 
malaria elimination strategies [23,62]. Finally, the poten- 
tially disastrous consequences of the rise and spread of 
artemisinin resistance requires that detailed and effective 
planning be implemented in preparation for containing 
and stemming any spread [56]. A basic assessment here 
were provided of prevalence-scaled travel from the four 
Southeast Asian countries where resistance has previ- 
ously been observed, but significant refinements of these 
estimates and modelling methods should be undertaken. 
These may include improved tracking and mapping of 
observed resistance and human movement patterns in 
Southeast Asia, as is being undertaken by the TRAC 
project [63], as well as scenario modelling of the risks of 
resistance escape to Africa or Latin America. Further, 



the incorporation of accessibility [64,65] and travel data 
[51,59] with drug use data (e g, [66]), prevalence infor- 
mation [40,41] and models [67], all undertaken within a 
probabilistic modelling framework (e. g., [6,60]), could 
aid in estimation of spread routes should resistance arise 
elsewhere. 

Additional files 



Additional file 1: The travel time/distance mask to extract the 
estimated typical maximum prevalence of Plasmodium falciparum/ 
Plasmodium vivax at the origin of travellers. Inset map: travel time to 
the nearest major settlement (population size >50,000). The global map 
of accessibility is obtained at [44]). Main map: each dot shows an 
airport location with a 50-km buffer around it, and the colours show 
the global P. falciparum prevalence map [40] masked by the global 
travel time map with a threshold value of less than two hours. These 
two-hour and 50-km thresholds were used to assign prevalence values 
to airports (see main text). 

Additional file 2: Air travel network communities weighted by 
directed estimates of passenger flow. Airports with the same 
community membership (indicated by the same colour) display stronger 
links in terms of likely movement volume between them than to airports 
in other communities. The movement volume is extracted from Huang et 
al's [36] modelled passenger flow matrix. 

Additional file 3: Communities for all possible connections 
originating from Plasmodium falciparum/Plasmodium v/Vax-endemic 
areas. A) P. falciparum multilevel membership; B) P. vivax multilevel 
membership. These two maps show directly connected, one-transfer and 
two-transfer airports from endemic areas. The inset maps present close- 
up views of the USA and western Europe. Airports with the same 
community membership (indicated by the same colour) display stronger 
links in terms of likely movements of infections between them than to 
airports in other communities. Note that in the P. vivax map, two 
communities with less than ten airports are not shown. 

Additional file 4: Top ten airports based on estimated relative 
malaria importation and exportation rates. P. falciparum/P. vivax flow 
measures are calculated based on the incoming and outgoing numbers 
of passengers travelling internationally, scaled by the malaria prevalence 
at the origin of the routes in the case of importation, and at the airport 
listed in the case of exportation. The flows represent a relative measure 
of infection movement and are not designed to represent actual number 
of infections. 

Additional file 5: Spatial distributions of airports with Plasmodium 
falciparum/Plasmodium vivax betweenness centrality scores. A) Airports 
with normalized betweenness scores >0 from P. falciparum-endemlc areas, 
weighted by the P. fa/c/pari/m-prevalence weighted passenger flows. B) 
Airports with normalized betweenness scores from P. wVox-endemic areas, 
weighted by the P. vivax prevalence-weighted passenger flows. Details on 
the betweenness metric are provided in the main manuscript. 

Additional file 6: Spatial distributions of airport nodes in 
elimination countries [55] weighted by incoming international 
Plasmodium falciparum/Plasmodium vivax flows. A) Airports in 
countries with elimination objectives shown with dot size scaled to 
match total incoming passenger flow weighted by the P. falciparum 
prevalence at the traveller origins. B) Airports in countries with 
elimination objectives shown with dot size scaled to match total 
incoming passenger flow weighted by the P. vivax prevalence at the 
traveller origins. These two figures highlight the relative risks of infection 
importations through air travel for each country with malaria elimination 
objectives [55], 

Additional file 7: Top ten Plasmodium falciparum routes and 
Plasmodium vivax routes to non-Asian destinations from the Great 
Mekong subregion, as defined by malaria-prevalence scaled 
passenger numbers. Values with (*) are returned as the adjusted largest 
flows between two connection flights. 
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